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Correlation measurements indicate that excess two particle correlations extend over causally disconnected ra- 
pidity ranges. Although, this enhancement is broad in relative rapidity r) = rji — r/2, it is focused in a narrow 
region in relative azimuthal angle <j> = <j>\ — (j>2. The resulting structure looks like a ridge centered at r\ = <j> = 0. 
Similar ridge structures are observed in correlations of particles associated with a jet trigger (the hard ridge) and 
in correlations without a trigger (the soft ridge). The long range rapidity behavior requires that the correlation 
originates in the earliest stage of the collision, and probes properties of the production mechanism. Glasma 
initial conditions as predicted by the theory of Color Glass Condensate and provide a and early stage correlation 
that naturally extends far in rapidity. We have previously shown that the soft ridge is a consequence of particles 
forming from an initial Glasma phase that experience a later stage transverse flow. We extend this work to 
study the ridge dependance on the pt of the correlated pairs. We then determine the soft contribution to the 
hard ridge. 



1. Introduction 



2. Glasma Correlations 



Correlation measurements at RHIC show an en- 
hancement as a function of relative rapidity r/ = 
Vi ~~ an d relative azimuthal angle (f> — 4>i — 4>2- 
Centered at r] = <fi = 0, the ridge gains its title be- 
cause the measured structure is broad in 77 and narrow 
in <p. Both triggered and un-triggered measurements 
show a ridge like structure, but have been considered 
separate phenomena. The hard ridge measures the 
correlated yield of associated particles per jet trigger 
PQ, and the soft ridge measures the number of cor- 
related particle pairs per particle [2]. We propose a 
common explanation for these phenomena based on 
particle production in an early Glasma stage followed 
by radial flow [3]. 

In [3] we extend the model of Ref.|3] to incorpo- 
rate jet production and address the soft and hard 
ridges. Long range rapidity correlations and the in- 
sight they provide on early time dynamics are our 
driving concerns [4J |5]- We take particle production 
to occur through a Glasma state, and emphasize how 
computed Glasma correlations can affect ridge mea- 
surements. In Ref. jlj we found excellent agreement 
with the peak amplitude and azimuthal width shown 
in current Au+Au data. In the next section we ex- 
tend this work here to include Cu+Cu systems. We 
then extend the model of Ref. [I] to address varying p t 
ranges so that we may address the hard ridge. In Sec. 
[3] we add a contribution of jets following the model of 
Ref. [5J. We extend that model to compute both the 
strength and azimuthal dependence of the jet contri- 
bution to the hard and soft ridges. We then combine 
the flow and jet effects and find that correlations of 
thermally produced pairs constitute a significant con- 
tribution to the triggered measurement. We then dis- 
cuss how experiments might distinguish the different 
contributions. 



The theory of Color Glass Condensate (CGC) pre- 
dicts an early Glasma stage in a high energy colli- 
sion in which particles are produced by strong longi- 
tudinal color fields. As nuclei collide, the transverse 
fields of each nucleus are instantaneously transformed 
into longitudinal fields that are approximately uni- 
form in rapidity. These fields are essentially random 
over transverse distances r t larger than the saturation 
scale Q7 1 , where Q s ~ 1 — 2 GeV. We can think of 
such field configurations as consisting of a collection 
of longitudinal flux tubes. Flux tubes are ubiquitous 
in QCD-based descriptions of high energy collisions. 
In the Glasma they are closely packed and not strictly 
distinct due to saturation. They are, however, uncor- 
rected for r t > Q7 1 . This is their essential feature for 
this work. The Glasma changes to plasma as particles 
form from the fields and thermalize. 

In the saturation regime, the number of gluons in a 
rapidity interval Ay is 



N = {dN/dy)Ay ~ a s 



Q S R A , 



(1) 



where Ra is the nuclear radius and a s is the strong 
coupling constant at the saturation scale Q s [7]. We 
understand as the number of flux tubes K ~ 
(Q s Ra) 2 times the density of gluons per flux tube 



The scale of correlations is set by 



K = 



(N 2 ) - (N) 2 - (N) 

W 2 • 



(2) 



where the brackets denote an average over collision 
events. This quantity vanishes for uncorrelated glu- 
ons, for which multiplicity fluctuations are necessarily 
Poissonian. In Ref. [I] we argued that the Glasma 
correlation strength is 1Z oc (K)^ 1 = (Q S R)~ 2 , and 
found that the Glasma contribution to correlations is 



KdN/dy ~ a s (Q s y 



(3) 



2 
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a result consistent with calculations of Dumitru et al. 
in Ref. |5j. Equations ([!]) and |3| constitute initial 
conditions for the hydrodynamic evolution of the sys- 
tem. 

These Glasma initial conditions affect the final state 
correlations in several ways. First, particles emit- 
ted from the same tube share a common origin that 
is localized to a very small transverse area, since 
Q s <C Ra- Second, the flux tubes correlate particles 
over a large pseudorapidity range. Some of the flux 
tubes can stretch across the full longitudinal extent of 
the system at times < 1 fm. These flux tubes rapidly 
fragment. At later times, the particles they produce 
can be separated by large longitudinal distances de- 
pending on their momenta. Consequently, subsequent 
scattering and hydrodynamic evolution cannot erase 
their correlations - they are causally disconnected. 
Third, the strength of the correlation depends on the 
number of flux tubes. The number of tubes depends 
on the centrality and energy of the collision, as well 
as the transverse area of the tube, all of which, in 
turn, depend on Q s . Finally, as a result of the com- 
mon origin, particles coming from the same tube must 
have the same initial radial position and feel the same 
effects from flow, independent if their rapidity. 

During the Glasma phase, flux tubes thermalize 
into partons and pressure builds as the systems moves 
toward an equilibrated Quark Gluon Plasma. Par- 
tons initially localized in tubes are now localized in 
small fluid cells with a uniform azimuthal distribu- 
tion of particles. Following a Hubble-like expansion 
of the system, the transverse fluid velocity takes the 
form 7 t v t = Xr t . All of partons in a fluid cell are 
boosted radially depending on their initial radial po- 
sition. Consequentially, partons in any given fluid cell 
gain transverse momentum in the radial direction and 
the relative angle between any two momentum vectors 
in that cell becomes smaller. Furthermore, fluid cells 
at a larger radial position have a larger final trans- 
verse velocity and the relative angle between parton 
momentum vectors is narrower. This angular narrow- 
ing depends only on the initial radial position and the 
small transverse area of the flux tube source. In our 
simplistic view, all flux tubes are uniform in rapid- 
ity and extend to the same longitudinal length. This 
provides for a longitudinally uniform fireball that ex- 
periences the same radial flow at every longitudinal 
position. In this way, at every rapidity, angular corre- 
lations are enhanced in the same way and the initial 
state spatial correlations are both preserved through 
freeze out and represented in momentum correlations. 

Following [4] we define the momentum space corre- 
lation function at freeze out as 

KP1-P2) = P2(Pl,P2) - Pl(Pl)Pl(P2) (4) 

where p 2 (pi,p 2 ) = dN/dyid 2 ptidy 2 d 2 p t2 is the pair 
distribution, and pi(p) = dN/dyd 2 p t is the single par- 
ticle spectrum. We describe the effect of flow using 



the familiar blast- wave model [H EE EU EEJ EH- In 
this model, the single particle spectrum is pi(p) = 
/ /(x,p)dr, where /(x,p) = (2^)" 3 exp{-p"^/T} 
is the Boltzmann phase-space density and dT = p^da^ 
is the differential element of the Cooper-Frye freeze 
out surface. We assume a constant proper time freeze 
out, so that dT = Tpm t cosh(y — r))dr)d 2 r t , where 
-q = (1 /2) ln((< + z)/(t — z)) is the spatial rapidity. 

We argue in Ref. jlj that the final-state momentum 
space correlation function is 

1 \ f 1 n/( x i.Pi) /( x 2,P2) rr jr, 

r-(pi,p 2 ) = / c ( x i: x 2) j — v , — ^dT 1 dT 2 , 

J ni(xi) ni(x 2 ) 

(5) 

where ni(x) = J f(x,p)dT. The spatial correlation 
function c(xi,X2) depends on the Glasma conditions 
as follows 



c(xi,x 2 ) = K6(r t )p FT (R t ), 



(6) 



where r t = r tl — r t2 is the relative transverse position, 
and Rt = (r t i + r t2 )/2 is the average position. The 
delta function accounts for the fact that Glasma cor- 
relations are highly localized to r t < Qj 1 . The factor 
p FT (R t ) describes the transverse distribution of the 
flux tubes in the collision volume, which we assume 
follows the thickness function of the colliding nuclei. 

We comment that the form of ^ holds as long as 1Z 
is unmodified from its initial Glasma value by particle 
production and hydrodynamic evolution. This is only 
strictly true as long as a) subsequent evolution doesn't 
change the relative number of particles in the rapid- 
ity interval of interest; and b) the number of observed 
hadrons is proportional to the initial number of glu- 
ons. Causality prevents these effects from altering 1Z 
for truly long range correlations, \r)i — ?7 2 | > 1 — 2. This 
would hold for smaller rapidities in Glasma theory as 
long as boost invariance is a reasonable approxima- 
tion, since dN/dy is then a hydrodynamic constant 
of motion in each event. Moreover, assumption (b) is 
common in Glasma/CGC calculations. On the other 
hand, for 1 771 — 772 1 < 1 — 2, the experimental dN/dy 
is not flat and will change with time due to particle 
diffusion and number changing processes; see |13j . For 
now, we will assume that 1Z is constant and defer the 
hydrodynamic modification for later work. 

The analysis in Ref. |4J focused on a measurement of 
the near side peak of the soft ridge in 200 GeV Au+Au 
using the observable Ap/^/p = (p sib - p ref )/^fp^ 
[2J. The quantity p S ib(4>,ri) represents the distribu- 
tion of "sibling pairs" from the same event, as a 
function of relative pseudorapidity and relative az- 
imuthal angle, and is comparable to our p 2 . The 
quantity p re f((j),r]) represents uncorrelated pairs from 
mixed events and is equivalent to the square of our 
p\. The difference (p S ib — Pref) is a measure of cor- 
related pairs and is comparable to the integral of Q 
over the transverse momenta and average azimuthal 
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angle <& = (</ 
the quantity 



<\>2)j2. It is convenient to compute 



A/3 _ / r(p ll p 2 )ptiPt2dptidp t 2d$ 

Pref J Pl(Ptl)pi(Pt2)d 2 ptld 2 p t 2 



(7) 



which is independent of the overall scale of the mul- 
tiplicity. We emphasize that this quantity includes 
correlated pairs in which both particles can have any 
momentum, while measurements of the hard ridge cor- 
relate particles from different pt ranges. We will ex- 
tend our approach to address such quantities below. 

To construct the observed quantity Ap/ ^/p re f, we 
notice that r(pi,p2) computed from eq. (ml is pro- 
portional to the correlations strength TZ pIT We can 
write 



Ap_ 

Pref 



HF(<j>), 



(8) 



where F(4>) is normalized such that J Q 27r F((j))d(f) = 1. 
The distribution F(<f>) depends only on the blast-wave 
parameters ^ym/T and v s , and represents the angular 
correlations of particles from flux tubes after hydrody- 
namic expansion. The factor TZ scales the strength of 
the correlations with both energy and centrality and 
determines the rapidity dependence (which is flat in 
this case). 

We now combine ^ and ([!]) to obtain the observed 
quantity 



Ap/Vp^7 = nTZdN/dy F(4>), 



(9) 



where we equate the factor TZdN/dy with ([3]), which 
accounts for all of the Glasma energy and centrality 
dependence. The scale constant k is independent of 
energy. As described in jl] we set k only for Au+Au 
200GeV collisions such that F(<j>) for the most cen- 
tral collisions is aligned with the most central data 
point. Although blast- wave parameters have some 
energy dependence, the Glasma factor ^ allows for 
strong agreement with the 62 GeV data without fur- 
ther adjustment of k. 

We now apply ^ and ^ to address new and forth- 
coming data for Ap/ y/p re f and the azimuthal width 
Cu+Cu for both 200 and 62 GeV as a function of 
centrality. No parameters are adjusted from pf - in 
that sense, these results are predictions. The Q 2 S has 
a dependence on the density of participants p par t that 
is determined in Ref. [TJ. The saturation scale in 
central Cu+Cu is then 



s (Cu) 5 



= (Au) 



2 P P art(Cu, central) 
Ppart (Au, central) 



(10) 



After this scaling, we take the relative centrality de- 
pendence of Q s to be the same in the Cu and Au sys- 
tems. This assumption can be refined by measuring 
the centrality dependence of dN/dy in Cu+Cu. Sim- 
ilarly, we obtain the blast wave parameters T and v s 
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Figure 1: Cu+Cu amplitude as a function of centrality 
(y = 2N bln /N part ) for both 200 and 62 GeV. Dashed lines 
represent the blast wave blast-wave integrations F((j>) nor- 
malized to 200 GeV Au+Au as in [4] . Solid lines represent 
the inclusion of the CGC scaling, TZdN/dy, following |9| 
with the inclusion of ( 10 1 . The parameters are unchanged 
from Ref. g] 



in Cu+Cu from the Au+Au values by assuming that 
they scale with the number of participants. These 
assumptions can be refined as single particle spectra 
measurements as a function of centrality become avail- 
able. 

Figure [T] shows our prediction for the soft ridge am- 
plitude in Cu+Cu 200 and 62 GeV systems as a func- 
tion of centrality. In both panels, the dashed line rep- 
resents the blast-wave amplitudes F((f> = 0) — F((f> — 
7r) normalized to fit Au+Au at 200 GeV in Ref. 0]. 
The dashed lines are included to show how the en- 
ergy and system-size dependence affects the calcula- 
tion when the Glasma dependence ^ is omitted. The 
solid lines are the result of including the Glasma scal- 
ing ([3]) adjusted by (JToJ) . The error band represents a 
10% uncertainty in the blast- wave parameters plus an 
additional uncertainty in the parameterization of Q s 
that increases with decreasing centrality. 

In Fig. [2] we show the soft ridge azimuthal width 
in Cu+Cu systems compared to previously published 
Au+Au result [I]. We have also included preliminary 
STAR measurement of the soft ridge in Au+Au width 
[2] . In [4] we find that the azimuthal width of the near 
side peak of the soft ridge in Au+Au is due to radial 
flow, is constant with a change in energy, and is rel- 
atively uniform with change in centrality. The error 
band is representative of the uncertainty of fitting an 
offset gaussian to the angular calculation. Since the 
azimuthal width is completely determined by radial 
flow, which depends completely on the choice of cen- 
trality and blast-wave parameterizations, and all of 
those parameters have remained unchanged, we cal- 
culate the same enhancement in the width for Cu+Cu 
as Au+Au. The black line extending to v = 6 is the 
Au+Au result from [4], and the overlaid blue line ex- 
tending to v = 4 with the hatched error band is the 
Cu+Cu result. Again, as with the Au+Au result, the 
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Figure 2: Comparison of the previously published angular 
width calculations for Au+Au (black line) with the angu- 
lar widths for Cu+Cu systems (blue line, hatched error 
band) compared with preliminary STAR Au+Au 200 and 
62 GeV data. Width calculations remain independent of 
energy and nearly independent of system 



Cu+Cu result is independent of energy since the mea- 
sured transverse expansion does not depend on energy. 

A key feature of the flow-based descriptions of the 
ridge is that it is the angular width, oy, of correlated 
pairs decreases as the mean p t of the pair increases. 
The greater the radial boost given to a fluid cell, the 
narrower the relative angle between the momentum 
vectors of particles in that cell. A very high momen- 
tum correlated pair is more likely to have come from 
a fluid cell that received a very large transverse boost . 
To study whether this effect is present, we compute 
Ap/y/Pref for pairs of p t >Pt, m in- As pt, m i n increases, 
we also expect the amplitude of the ridge to decrease, 
since it is more difficult to find higher p t bulk particles. 
We therefore compute 



Ap 

Pref 



Pt2,max Ptl,,nax 

I J r(pi,p 2 ) 

Pt2,min Ptl,min 



Ptl, moi 



Pt1,r, 



Ptl.min Pt2,min 



TlF(cj);pti, mm ,pti,, 



,Pt2,. 



,Pt2, 



(11) 



where the integration measures are the same as in ^ . 
To obtain the measured ratio we write 



. oo 



poo 

) Jptl, 



Pi 



Pt,min, OO-iPt,' 
Ap/^/P^ef F{(j)) J °° Pl 

(12) 

For increasing values of pt, m in we calculate the cor- 
relation function as before using ( [Tl) and Q and 
find that the azimuthal width does indeed decrease as 
shown in the upper panel of Fig. [3j 

The upper panel of Fig.[4]shows the correlation am- 
plitude vs. centrality for different choices of pt,min- 
The amplitude decreases with increasing pt, m in be- 
cause the number of particles contributing to correla- 
tions is reduced. The lower panel of Fig. [4] shows the 
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Figure 3: Top panel shows the angular width for Au+Au 
systems with increasing minimum pt,min limits compared 
with preliminary STAR Au+Au 200 GeV data. The 
dashed line represents the pt,min=0 calculation shown in 
Fig. [2] The lower panel shows the azimuthal width for 
most central collisions vs the pt,min limit for both Au+Au 
and Cu+Cu 200 GeV. 



amplitude of Ap/ \jp r ef for the most central collision 
as a function of the choice of pt, m in for Au+Au and 
Cu+Cu at 200 GeV. Similarly, the blue curve in the 
lower panel of Fig. [3] represents the azimuthal width of 
the soft ridge in most central collisions as a function 

of Pt,min • 

We see that the azimuthal width of the hard ridge 
is smaller than that of of the soft ridge, but as the 
Pt,min limit of the soft ridge is increased, the ampli- 
tude of the correlations drops and the azimuthal width 
narrows. In the pt range of the hard ridge, it appears 
that the azimuthal width could be narrow enough, but 
the amplitude is not directly comparable. To under- 
stand this difference, we must understand the differ- 
ences in the two measurements. The most significant 
difference is the choice of the momentum range of the 
correlated particles. The hard ridge measurement an- 
alyzes the yield of associated particles per jet trigger 
where the associated particle pt range and the trigger 
range do not overlap. The soft ridge measurement, 
however, finds the number of correlated pairs per par- 
ticle where both particles are in the same range with 
p t above minimum bias. The normalization of the soft 
ridge is found by taking the square root of the uncor- 
related pair reference spectrum. 

STAR measures the hard ridge, or yield of associ- 
ated particles per jet trigger, for Au+Au 200 GeV for 
3 < Pt.trigg < 4 GeV with 2 < p t ,assoc < 3 [T]. Iden- 
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Figure 4: Top panel: Au+Au 200 GeV amplitude cal- 
culations for increasing pt,min limits. The dashed line is 
calculation with pt,min = 0. Lower panel: the soft ridge 
amplitude for most central collisions plotted as function of 
the pt,min limit for both Au+Au and Cu+Cu 200 GeV. 



tifying p t \ with the trigger range and the associated 
range with p t 2, we calculate Ap/y/p re f and transform 
to yield by 



Yield = 



( 



J PiiPti) 



(13) 



At higher ranges of pti 2 the contribution from jets 
should become more significant. It is important 
therefore to know the relative contribution of ther- 
mal particles and jet particles. As will be discussed 
in more detail later, we decompose the total par- 
ticle spectrum into thermal bulk and jet fractions, 
and to obtain the contribution of bulk correlations 



to the hard ridge, we multiply ( 13 ) by the bulk frac- 
tion J piipti)/ j Ptot(pti) where Jp t ot(pti) is the total 
number of particles in the range of pti- 

The blue curve in Fig. [5] represents the contribu- 
tion to the hard ridge from only thermal bulk pairs. 
As can be seen on the figure, bulk-bulk correlations 
contribute significantly to the amplitude of the trig- 
gered measurement, but seems to have a somewhat 
narrow profile in azimuth. It was shown in [fjj that a 
jet acquires angular correlations with flowing matter 
due to quenching, but the width of the correlation is 
wider than the data. The contribution of jet corre- 
lations with bulk particles could make up the differ- 
ence between the blue curve in Fig. [5] and the data by 



increasing both the amplitude and the width of the 
calculation. 

In the next section we combine a theory angular 
correlations from 6 with spatial correlations of jets 
and flux tubes to obtain a jet-bulk contribution to the 
hard ridge. 



3. Jets, Glasma, and Correlations 

As the pt of correlated particles is increased, the 
contributions from jets should become prevalent, par- 
ticularly for small rj. At small rapidity differences, 
correlations of jet particles with fragments should be 
large, but restricted to the size of the jet cone, and 
the transfer of momentum from jet particles to bulk 
particles is causally limited to ~ 1 — 2 units in rapid- 
ity. The existence of correlations with jet particles at 
larger r\ would require a correlation early in the colli- 
sion that remains through the longitudinal expansion 
of the system and is still present at freeze out. Both 
the hard collisions and flux tubes are made in the ini- 
tial moments of the nuclear collision. Assuming that 
the entire overlap region of the colliding nuclei is in 
the saturation regime, flux tubes would fill the colli- 
sion volume and a jet formed at any transverse posi- 
tion would be accompanied by a flux tube at the same 
position. Since the flux tube extends to large rapidi- 
ties, the correlation of particles from the tube and the 
jet can extend to large n. Angular correlations arise 
since particles from the tube acquire a radial trajec- 
tory from flow as before, but the jet trajectory has a 
bias in the radial direction due to quenching [6 . 

We construct a distribution of jets as follows. We 
assume that jets are produced with a hard scattering 
rate /o(Pi) that is independent of position, multiplied 
by a spatial profile P ro d oc (1 — r\j R\) that is roughly 
proportional to the density of binary collisions. The 
phase space density of jet particles is then 

/./(xi,pi) = /o(pi)-Pprcd(n).S(r-i,0i), (14) 

where S is the survival probability of a jet due to jet 
quenching. In practice /o(pi) cancels in Ap/p re f so 
we need not specify it. We follow Ref. [H] and take 



S(ri,<f>i) = exp(-L(ri,<t>i)/l a bs), 



(15) 



where l a i,s — 0.25/rn is the jet attenuation length. 
The survival of the jet depends on the path it takes 
out of the medium 



R 2 



r?sin 2 ( 



i) - n cos(r/>i). (16) 



The path ( 16 ) is the distance a jet would have to travel 
out of a circular transverse area at an angle <f)\ with 
respect to the radial vector pointing to its position 
of production r± [5]. In view of (16 1 and (15), the 



shortest path, which is the path a jet is most likely to 
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Figure 5: Angular profile of the jet triggered ridge in a ra- 
pidity range away from the jet peak. The solid black line 
combines of long range correlations. Bulk-bulk (blue line) 
and jet-bulk (dashed red line) contributions are shown sep- 
arately. The jet fraction is determined by p 3 = 1.25 GeV; 
bulk-bulk correlations make up ~ 75% of the total ampli- 
tude. 



survive, is one that is radially outward from its posi- 
tion of production. Although a jet production is at a 
minimum on the surface, r w Ra, with little material 
for the jet to pass through, the survival probability 
is maximum in any direction (not pointing into the 
volume). The resulting angular correlations are weak, 
since only a small fraction of phase space is restricted 
by quenching. The largest probability of production 
would occur at the center where L ss Ra in all direc- 
tions, but there are no correlations in this case since 
quenching would be maximum. More concisely, jets 
are less likely to be produced at the surface and have 
a wide angular distribution. The most jets are pro- 
duced the center and would have the narrowest cor- 
relation with radially flowing particles, but have the 
highest probability of being quenched. As the pro- 
duction point of the jet moves from the center toward 
the surface, the probability of production decreases, 
the probability of survival increases, and the angular 
correlation with radially flowing particles widens. In- 
tegration over all possibilities determines the width of 
the correlations. 

The calculation follows the analysis of Q , but with 
the first particle from a jet and the second from a flux 
tube, so that 



! (Pl,P2) = / c JB (xi,x 2 



, /■/(*!, Pi) /b(*2,P2) 

«ij(xi) ni B (x 2 ) 



dT 1 dT 2 
(17) 

The correlation function c JB (xi, x 2 ) requires that a jet 
and a bulk particle must come from the same radial 
position. The rest of the equation accounts for the 
different the spectra for jet and bulk particles. We 



»(xi,x 2 ) 



K JB (Nj) 

n(N B ) 



c(xi,x 2 ), 



(18) 



where c is given by ([6]). 

To relate the correlation strength 1Z JB to the bulk 
correlation strength 1Z discussed earlier, we assume 
that the hard scattering rate is independent of the 
flux tube dynamics. Recall that the bulk quantity 
7Z in (2) is related to the number of flux tubes; see 
([3]). If we take the fraction of jet and bulk particles 
per flux tube to be independent of the number of flux 
tubes, we can write (Nj) = a(N), where a and (3 
may depend on momentum, but do not vary event by 
event. We then follow [H] to find (N B ) = f3{N), and 
(NjN B ) = af3(N(N — 1)), so that 



n 



(Nj){N B ) 
af3(N(N-l)) -a{N)l3{N) 



a(N)P(N) 
(N(N-l)) - (AO 2 

w 



K. 



(19) 



Therefore, the addition of jets to the total multiplicity 
doesn't change the correlation strength. In essence, 
the beam jet associated with the hard process is just 
another flux tube in the high density Glasma state. 
We can now rewrite (J9J) for jet-bulk correlations as 

^Pju/Vf^f = KKdN jet /dy F JB (0). (20) 

Our calculation of yield and implementation of lower 
Pt limits follows (13 1 and (ITT]) but, this time, we scale 
by the jet fraction/ p\,j{pn}/ J pi,tot(pti)- 

In order to compute the amplitude and azimuthal 
width of the hard ridge in Fig. [5] we must determine 
the relative contributions from both bulk-bulk and jet- 
bulk correlations. The measured pt spectrum follows 
an exponential behavior at low p t and a power law be- 
havior where jets play a larger role, see e.g. Ref. [T5] , 
Our blast wave formulation describes the exponential 
behavior of the low p t spectrum well. The scale p s at 
which the spectrum begins to deviate from exponen- 
tial behavior is proportional to Q s in Glasma theory, 
but the proportionality constant is not known. This 
introduces a free parameter - p s at y/s = 200 GeV - 
that we fix below. We then find the number of jet par- 
ticles by taking the difference between the total num- 
ber of particles and the number of thermal particles 
Pi,j — P\,tot — Pi.b- We take pi,tot from the measured 
spectrum in Ref. [15 and p\ t B from the blast wave 
calculation with the appropriate normalization. 

We now calculate the combined effect of Glasma, 
flow, and jets on the correlation function. Adding the 
bulk-bulk and jet-bulk contributions, we obtain 



= ^F B 

y/Pref dy 



I piAPt) 

f Pl.tot(Pt) 
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dN 3et / n S Pi APt) 



+K n^F JB {4>) 



J Pl,tot(Pt)' 



(21) 



In order to compare to the yield of associated parti- 



cles in the hard ridge, we combine ( |21| and (13 1 with 
the appropriate integration limits. We find that the 
agreement with the data in Fig. [5] requires p s — 1.25 
GeV. The dashed red curve in Fig. [5] represents the 
contribution to the yield from correlations of jet and 
bulk thermal particles. This contribution is too wide. 
On the other hand, the bulk-bulk correlation function 
describing the effect of flow alone, which is given by 
the blue curve, is too narrow and the computed peak 
height is too small. The combination of the two effects 
shown as the black curve in Fig.|5]gives nice agreement 
with both the amplitude and azimuthal width. 

We now compare hard and soft ridge measurements 
directly by computing the momentum dependent cor- 
relation function (21 1. At low pt, m in the contribution 
from jets is negligible, therefore the amplitude and 
width of (21 1 is determined by the bulk-bulk term. As 
the Pt, m in is increased, both the amplitude and the az- 
imuthal width of the bulk-bulk term decreases, while 
the amplitude of the jet-bulk term increases. The az- 
imuthal width of jet-bulk correlations is roughly inde- 
pendent of pt,min', the growth is due to the growth of 
the jet fraction. Jet-bulk correlations should become 
a more significant fraction of the total as Pt, m in is m ~ 
creased, and the azimuthal width of the ridge should 
increase toward the jet-bulk width. 

We emphasize that the decrease of the bulk-bulk 
contribution to the amplitude of Ap/^/p and o~ r with 
increasing p t in Figs. [3] and [4] is a direct consequence 
of transverse flow and, consequently, is a firm pre- 
diction if the jet-bulk contribution is neglected. The 
role of jets and other phenomena like recombination 
are less clear. We have chosen the model of Ref. [B] 
because it relies only of the well-studied phenomena 
of jet quenching. Our calculations using this model 
predict that the width would increase for higher p t 
ranges. 



4. Summary 

We have shown that long range Glasma correla- 
tions contribute significantly to both the hard and 
soft ridges. Bulk-bulk correlations arise due to parti- 
cle production from a Glasma state followed by trans- 
verse flow [5], OH]. The Glasma initial conditions al- 
low us to extend to study of Cu+Cu systems at both 
200 and 62 GeV. We see in Fig [I] and Figj2j that we 
maintain good agreement with current data. To study 
the soft contribution to the hard ridge, we study the 
momentum dependence of bulk-bulk correlations and 
find that they are still significant, even in the p t ranges 
of the hard ridge measurement. To explore the effect 



of jet production and quenching on the ridge we in- 
cluded a model of jet-bulk correlations following [6 . 
As seen in Fig [5] we find that we cannot explain the 
magnitude and width of the hard ridge without a sub- 
stantial soft component. 
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